Search for: All records

Abstract Historical accounts suggest that Euro-American agricultural practices (post–1850 CE) accelerated soil erosion in the Paleozoic Plateau of the Upper Mississippi River Valley (USA). However, the magnitude of this change compared to longer-term Late Pleistocene rates is poorly constrained. Such context is necessary to assess how erosion rates under natural, high-magnitude climate and eco-geomorphic change compare against Euro-American agricultural erosion rates. We pair cosmogenic 10Be analyses and optically stimulated luminescence (OSL) ages from samples of alluvium to build a paleoerosion-rate chronology for Trout Creek in southeastern Minnesota (USA). Erosion rates and their associated integration periods are 0.069–0.073 mm yr−1 (32–20 ka), 0.049 mm yr−1 (28–14 ka), and 0.053 mm yr−1 (14–0 ka). Based on previous studies, we relate these rates to (1) the transition from forest to permafrost at the onset of the Last Glacial Maximum, (2) the decline of permafrost coupled with limited vegetation, and (3) climate warming and vegetation re-establishment. These pre-settlement erosion rates are 8× to 12× lower than Euro-American agricultural erosion rates previously quantified in the region. Despite a limited sample size, our observed rapid increase in erosion rates mirrors other sharply rising anthropogenic environmental impacts within the past several centuries. Our results demonstrate that agricultural erosion rates far exceed climate-induced erosion-rate magnitude and variability during the shift from the last glaciation into the Holocene. 

Abstract The Channeled Scabland of eastern Washington (USA) was formed by outburst floods from glacial Lake Missoula. Despite chronological advances, the timing of erosion in the main flood channels is unresolved. In particular, it is still uncertain whether upper Grand Coulee, the largest canyon in the Channeled Scabland, was incised during or prior to the last glaciation. We report 10Be exposure ages from erratics in upper Grand Coulee, glacial Lake Columbia, and surrounding flood routes. Flood-transported boulders on the high-elevation east rim of Grand Coulee date to ca. 17–15 ka. Ages from boulders on the floor of Grand Coulee indicate later flooding at ca. 14 ka, which post-dated canyon incision and occurred after inundation of the Telford-Crab Creek scabland at ca. 15–14.5 ka. Prior hydraulic modeling and dating suggest the entrance to Grand Coulee was blocked by rock and that canyon incision was incomplete at ca. 17 ka; hence, we interpret the 17–15 ka exposure ages on the east rim to coincide with flow over a retreating cataract during canyon incision. Our results indicate incision of Grand Coulee was completed between 17 ka and 14 ka. The short duration of canyon incision suggests that glacial Lake Missoula generated some of the most erosive outburst floods in Earth's history. 

Abstract Pleistocene outburst floods from the drainage of glacial Lake Missoula carved bedrock canyons into the Columbia Plateau in eastern Washington, USA, forming the Channeled Scabland. However, rates of bedrock incision by outburst floods are largely unconstrained, which hinders the ability to link flood hydrology with landscape evolution in the Channeled Scabland and other flood-carved landscapes. We used long profiles of hanging tributaries to reconstruct the pre-flood topography of the two largest Channeled Scabland canyons, upper Grand Coulee and Moses Coulee, and a smaller flood-eroded channel, Wilson Creek. The topographic reconstruction indicates floods eroded 67.8 km3, 14.5 km3, and 1.6 km3 of rock from upper Grand Coulee, Moses Coulee, and Wilson Creek, respectively, which corresponds to an average incision depth of 169 m, 56 m, and 10 m in each flood route. We simulated flood discharge over the reconstructed, pre-flood topography and found that high-water evidence was emplaced in each of these channels by flow discharges of 3.1 × 106 m3 s−1, 0.65–0.9 × 106 m3 s−1, and 0.65–0.9 × 106 m3 s−1, respectively. These discharges are a fraction of those predicted under the assumption that post-flood topography was filled to high-water marks for Grand and Moses Coulees. However, both methods yield similar results for Wilson Creek, where there was less erosion. Sediment transport rates based on these discharges imply that the largest canyons could have formed in only about six or fewer floods, based on the time required to transport the eroded rock from each canyon, with associated rates of knickpoint propagation on the order of several km per day. Overall, our results indicate that a small number of outburst floods, with discharges much lower than commonly assumed, can cause extensive erosion and canyon formation in fractured bedrock. 

Earth's drainage networks encode clues that can be used to decipher geologic and geomorphic history. Dendritic drainage patterns, the most common, typically form on approximately homogeneous bedrock. Variations in rock properties or lithologic structure can disrupt dendritic patterns and form, e.g., trellis or rectangular networks. Although textbooks include such lithological–drainage links, the mechanisms driving drainage reorganization via lithologic variability remain poorly understood. To cast light on this mystery, we study drainage patterns in post-glacial landscapes of the Upper Mississippi River Valley (UMRV). Pleistocene glaciers deposited till across parts of this region, burying a landscape of fluvially dissected sedimentary rock whose buried valley network differs from modern-day drainage patterns. As the current river network erodes and exhumes the bedrock, it comes to a geomorphic "decision point": Does it reorganize to recreate the paleodrainage network, or does it maintain its new drainage pattern? To understand this decision-making, we conducted idealized landscape evolution modeling experiments. Modeled landscapes that reintegrated more of the paleodrainage network exhibited higher tortuosity, measured by dividing the real flowpath length by shortest path-length to the outlet, and obtuse tributary-junction angles. We apply this metric to two adjacent landscape types in the UMRV: (1) never glaciated (Driftless Area, DA) and (2) formerly till-mantled (Driftless-style Area, DSA), and measure the basin-averaged tortuosity for sub-basins draining streams of order 1 through 7. Across the UMRV, tortuosity increases as the maximum stream order of the sub-basin increases. For each order, tortuosity is statistically higher in areas that had been previously buried and re-exhumed (DSA) than the DA, indicating that the rivers in the DSA have reintegrated the paleodrainage network since deglaciation. For the 1st and 2nd order sub-basins, the mean basin-averaged tortuosity in the DSA is ~1-2% higher than the DA (p-value < 0.01) and ~10-14% higher (p-value < 0.01) in the 6th and 7th order sub-basins. Our analysis suggests that a drainage-based metric, tortuosity, can identify landscapes where lithological heterogeneity or structure plays a dominant role in drainage organization. 

For landscapes to achieve a topographic steady state, they require steady tectonic uplift and climate, and a bedrock that is uniformly erodible in the vertical direction. Basic landscape evolution models predict that incising drainage networks will eventually reach a static geometric equilibrium – that is, the map-view channel pattern will remain constant. In contrast, natural rivers typically incise through heterogeneous bedrock, which can force reorganization of the drainage structure. To investigate how lithological variability can force landscape reorganization, we draw inspiration from formerly glaciated portions of the upper Mississippi Valley. In this region, depth-to-bedrock maps reveal buried dendritic river networks dissecting paleozoic sedimentary rock. During the Pleistocene, ice advance buried the bedrock topography with glacial till, resurfacing the landscape and resetting the landscape evolution clock. As newly formed drainage networks develop and incise into the till-covered surface, they exhume the buried bedrock topography. This then leads to a geomorphic "decision point": Will the rivers follow the course of the bedrock paleodrainage network? Or will they maintain their new pattern? Using a numerical landscape evolution model, we find that two parameters determine this decision: (1) the contrast between the rock erodibility of the glacial till (more erodible) and of the buried sedimentary rock (less erodible) and (2) the orientation of the surface drainage network with respect to the buried network. We find that as the erodibility contrast increases, the drainage pattern is more likely to reorganize to follow the buried bedrock valleys. Additionally, as the alignment of the two networks increases, the surface drainage network also tends to restructure itself to follow the paleodrainage network. However, when there is less contrast and/or alignment, the surface drainage pattern becomes superimposed on the bedrock topography, with streams cutting across buried bedrock ridges. Our results agree with field studies demonstrating that variability in erodibility exerts a first-order control on landscape evolution and morphology. Our findings can provide insight into how lithologic variation affects surface processes, drives drainage reorganization, and creates geopatterns. 

Chemical weathering influences many aspects of the Earth system, including biogeochemical cycling, climate, and ecosystem function. Physical erosion influences chemical weathering rates by setting the supply of fresh minerals to the critical zone. Vegetation also influences chemical weathering rates, both by physical processes that expose mineral surfaces and via production of acids that contribute to mineral dissolution. However, the role of vegetation in setting surface process rates in different landscapes is unclear. Here we use 10Be and geochemical mass balance to quantify soil production, physical erosion, and chemical weathering rates in a landscape where a migrating drainage divide separates catchments with an order-of magnitude contrast in erosion rates and where vegetation spans temperate rainforest, tussock grassland, and unvegetated alpine ecosystems in the western Southern Alps of New Zealand. Soil production, physical erosion, and chemical weathering rates are significantly higher on the rapidly eroding versus the slowly eroding side of the drainage divide. However, chemical weathering intensity does not vary significantly across the divide or as a function of vegetation type. Soil production rates are correlated with ridgetop curvature, and ridgetops are more convex on the rapidly eroding side of the divide, where soil mineral residence times are lowest. Hence our findings suggest fluvially-driven erosion rates control soil production and soil chemical weathering rates by influencing the relationship between hillslope topography and mineral residence times. In the western Southern Alps, soil production and chemical weathering rates are more strongly mediated by physical rock breakdown driven by landscape response to tectonics, than by vegetation. 

Abstract Erosion degrades soils and undermines agricultural productivity. For agriculture to be sustainable, soil erosion rates must be low enough to maintain fertile soil. Hence, quantifying both pre-agricultural and agricultural erosion rates is vital for determining whether farming practices are sustainable. However, there have been few measurements of pre-agricultural erosion rates in major farming areas where soils form from Pleistocene deposits. We quantified pre-agricultural erosion rates in the midwestern United States, one of the world's most productive agricultural regions. We sampled soil profiles from 14 native prairies and used in situ–produced 10Be and geochemical mass balance to calculate physical erosion rates. The median pre-agricultural erosion rate of 0.04 mm yr–1 is orders of magnitude lower than agricultural values previously measured in adjacent fields, as is a site-averaged diffusion coefficient (0.005 m2 yr–1) calculated from erosion rate and topographic curvature data. The long-term erosion rates are also one to four orders of magnitude lower than the assumed 1 mm yr–1 soil loss tolerance value assigned to these locations by the U.S. Department of Agriculture. Hence, quantifying long-term erosion rates using cosmogenic nuclides provides a means for more robustly defining rates of tolerable erosion and for developing management guidelines that promote soil sustainability. 

Significance The glacial Lake Missoula outburst floods are among the largest known floods on Earth. Dozens of these floods scoured the landscapes of eastern Washington during the last Ice Age, from 18 to 15.5 thousand years ago, forming what is known as the Channeled Scabland. We explored how changes in topography due to the solid Earth’s response to ice sheet loading and unloading influenced the history of megaflood routing over the Channeled Scabland. We found that deformation of Earth’s crust played an important role in directing the erosion of the Channeled Scabland. 

ABSTRACT New findings about old puzzles occasion rethinking of the Grand Coulee, greatest of the scabland channels. Those puzzles begin with antecedents of current upper Grand Coulee. By a recent interpretation, the upper coulee exploited the former high-level valley of a preflood trunk stream that had drained to the southwest beside and across Coulee anticline or monocline. In any case, a constriction and sharp bend in nearby Columbia valley steered Missoula floods this direction. Completion of upper Grand Coulee by megaflood erosion captured flood drainage that would otherwise have continued to enlarge Moses Coulee. Upstream in the Sanpoil valley, deposits and shorelines of last-glacial Lake Columbia varied with the lake’s Grand Coulee outlet while also recording scores of Missoula floods. The Sanpoil evidence implies that upper Grand Coulee had approached its present intake depth early the last glaciation at latest, or more simply during a prior glaciation. An upper part of the Sanpoil section provides varve counts between the last tens of Missoula floods in a stratigraphic sequence that may now be linked to flood rhythmites of southern Washington by a set-S tephra from Mount St. Helens. On the floor of upper Grand Coulee itself, recently found striated rock and lodgement till confirm the long-held view, which Bretz and Flint had shared, that cutting of upper Grand Coulee preceded its last-glacial occupation by the Okanogan ice lobe. A dozen or more late Missoula floods registered as sand and silt in the lee of Steamboat Rock. Some of this field evidence about upper Grand Coulee may conflict with results of recent two-dimensional simulations for a maximum Lake Missoula. In these simulations only a barrier high above the present coulee intake enables floods to approach high-water marks near Wenatchee that predate stable blockage of Columbia valley by the Okanogan lobe. Above the walls of upper Grand Coulee, scabland limits provide high-water targets for two-dimensional simulations of watery floods. The recent models sharpen focus on water sources, prior coulee incision, and coulee’s occupation by the Okanogan ice lobe. Field reappraisal continues downstream from Grand Coulee on Ephrata fan. There, some of the floods exiting lower Grand Coulee had bulked up with fine sediment from glacial Lake Columbia, upper coulee till, and a lower coulee lake that the fan itself impounded. Floods thus of debris-flow consistency carried outsize boulders previously thought transported by watery floods. Below Ephrata fan, a backflooded reach of Columbia valley received Grand Coulee outflow of small, late Missoula floods. These late floods can—by varve counts in post-S-ash deposits of Sanpoil valley—be clocked now as a decade or less apart. Still farther downstream, Columbia River gorge choked the largest Missoula floods, passing peak discharge only one-third to one-half that released by the breached Lake Missoula ice dam.